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A method of resummation of truncated perturbation series, related to diagonal Pade approximants 

but giving results independent of the renormalization scale, was developed more than ten years ago 

(w^ \ by us with a view of applying it in perturbative QCD. We now apply this method in analytic QCD 

Cn , models, i.e., models where the running coupling has no unphysical singularities, and we show that 

t^rj' the method has attractive features such as a rapid convergence. The method can be regarded as a 

r—{ ' generalization of the scale-setting methods of Stevenson, Grunberg, and Brodsky-Lepage-Mackenzie. 

The method involves the fixing of various scales and weight coefficients via an auxiliary construction 

of diagonal Pade approximant. In low-energy QCD observables, some of these scales become some- 

\^ ' times low at high order, which prevents the method from being effective in perturbative QCD where 

the coupling has unphysical singularities at low spacelike momenta. There are no such problems in 

analytic QCD. 
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I. INTRODUCTION 



^ ' Extending the applicability of QCD from high energies, where it can be consistently treated by perturbation 

^sj methods, down to the low-energy regime is one of the main tasks of theoretical hadronic physics. A simple-minded 

^^ utilization of perturbation series is clearly forbidden, not just by the sheer size of the expansion parameter (the 
^^ running coupling parameter a(Q^) = as{Q'^)/TT at low momentum transfer Q^ = —q^), but even more so by the 
^^ , existence of unphysical (Landau) singularities of the coupling parameter in the complex Q^ plane, the singularities 
t ' which are inferred from the renormalization group equation when the corresponding beta function is expressed in 
^^ ', terms of a truncated perturbation series. These singularities are unphysical because they do not reflect correctly the 
~^ analytic properties of spacelike observables 'D{Q^), properties based on the general principles of local quantum field 

theories [l', l2l|. Consequently, the most straightforward procedure for applying QCD to low-energy quantities consists 
in removing this unwanted nonanalyticity by some kind of analytization of the coupling parameter a((5") i— > ,4i((5"). 
The analytic coupling parameter Ai{Q'^) can differ significantly from the perturbative one a{Q^) only at low momenta 
IQ^I ^ 1 GeV^. Several constructions of such analytic QCD models, i.e., of Ai{Q^), have been made during the last 
C^ fifteen years - starting from the seminal papers of Shirkov et al. [3-Q. For reviews of various types of analytic QCD 
models see Refs. [6|-l9|. On the other hand, handling the physics of hadrons at low energies by simply utilizing an 
appropriately modified, "analytized," coupling parameter (together with its higher order analogs) within perturbative 
approaches is a very ambitious task, since it implicitly rests on the assumption that the low-Q^ behavior of AilQ"^) 
can be defined in a way that all nonperturbative effects are effectively included - at least for inclusive quantities. 
Of particular interest here is the behavior of Ai{Q'^) for Q^ -^ 0, and this question was the subject of intensive 
studies during last years, based either on analytic methods (Schwinger-Dyson equations 10], Banks-Zaks expansion 
[111 . uM ) o^ on numerical lattice approaches [l3| ■ They have finally led to the strong suspicion of "freezing" of the 
coupling parameter near Q^ = 0. If one wants to go a step further, however, and specify Ai{Q^) for the whole range 
IQ^I < Ql^g {Qas denotes the momentum transfer where asymptotic freedom should start to dominate) such that all 
non-perturbative effects get included, one clearly has to utilize as much as possible external information, both on the 
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side of empirical constraints and on the side of general physical principles such as causality, unitarity, analyticity, 
asymptotic freedom, operator product expansion, renormalization scale and scheme independence, etc. 

Within the present paper we focus mainly on the analytical structure and on the renormalization scale (RScl) inde- 
pendence of the resulting physical quantities. We apply, in various analytic QCD models, a global (i.e., nonpolynomial 
in the coupling) RScl invariant resummation/evaluation method which we developed in the context of perturbative 
QCD more than ten years ago [IJ, ll^, and we compare this evaluation method with other methods. In Sec. HIl we 
recapitulate the aforementioned RScl invariant resummation method for spacelike observables (in perturbative QCD). 
The presentation this time is somewhat less formal and, perhaps, more intuitive. In Sec. Illll we describe the minimal 
adjustments needed for the method to be used in analytic QCD models. In that Section we also argue why we should 
expect our resummation method to work significantly better in analytic QCD than in perturbative QCD. In Sec. IIVI 
we apply the method to the evaluation of the derivative of the massless (vector) current-current correlation function, 
i.e., the Adler function, both in perturbative QCD and in various motivated analytic QCD models. First, the evalu- 
ations are made for the leading-/3o part of the Adler function, where we know the exact result within each analytic 
QCD model, so this case is used as a test case for our resummation method to rather high values of the order index 
M . Subsequently, we apply our method to the truncated series of the Adler function, where only the first three full 
coeflficients (beyond the leading term) are known. In Sec.|V]we summarize the results and present conclusions. 

II. RECAPITULATION OF THE METHOD 

In this Section we present the resummation method developed in Refs. |14ill5l| in a somewhat simpler and, perhaps, 
more intuitive way. We consider a massless spacelike physical observable V{Q^) whose perturbation series in powers 
of the perturbative QCD (pQCD) coupling a{Q^) = as(Q^)/7r 






(1) 



2U2M] 



is known up to ^ a , such that we are faced with the truncated perturbation series 'D{Q ) ^ 

2M~1 



m')T = -iQ') + E d^ <Q"y^" 



(2) 



j=i 



Here we have chosen the renormalization scale (RScl) ji^ to be equal to the physical scale Q^ of the process (/i^ = Q^ 
For a general RScl /x^, the full and the truncated perturbation series read 



V{Q%t = a(A^2)+^d,(MVQ')«(A^'y+' 

2A/-1 



(3) 
(4) 



i=i 



This truncated series has a residual RScl dependence due to truncation. The /x -dependence of dj {fJ- /Q ) is dictated by 
the /i^ -independence of the full series 2?((5^)pt and the /^^-dependence of a(^^) given by the well known renormalization 
group equation 



da{fi'^) 
dln/i^ 



= - E ^^-2 a(M'r = -(3oai^iY (1 + cia(/x2) + c^aifi^ 



(5) 



J>2 



where the right-hand side is the beta function /3(a), and we denoted Cj = f3j/f3Q. In particular, we obtain (we denote 
throughout: dj{l) = dj and do = d^lp? /Q"^) = 1) 



rfi(/^VQ') = rfi+/3oln(MVQ') , 



d2(A*VQ') = rf2+E 



2! 



ik ^k / M 



,/3ol-'=l^)'^-. + /3iln(^ 



dz{^lVQ^) = ds+J2 






fc=l 



2*ln|i^)+|Aln^(|, 



-t-;32ln 



(6) 
(7) 

(8) 



etc. Note that a(/i^) and dj^j^/Q"^) are not only renormalization scale (RScl) dependent, but also renormalization 
scheme (RSch) dependent (as are also dj = dj{l)), i.e., they are functions of /x^, C2 = ^2/^0, C3 = Ps/Po, etc. The 
RSch dependence of dj{fi^/Q'^) and dj involves C2,...,Cj (when j > 2). The first two coefficients /Sq and /3i are 
universal in the mass independent schemes: /3o — (f 1 — 27i//3)/4, /3i — (f02 — 38ny/3)/16. 

In the following we will mainly be interested in the RScl dependence of the different (perturbation) series. Therefore, 
it will prove advantageous to use logarithmic derivatives of the pQCD coupling a instead of powers a". Specifically, 
we introduce^ 

_ (~i)" dMQ') (.. 

""+^(^^= (3^Snl d(lnQ2)n- (9) 

and reorganize the (truncated) perturbation series ©-Q into the "modified (truncated) perturbation series" (mpt) 



i?(0')mpt - a{^i') + Y,dJ{^lVQ')a,+,{^l^), (lo) 

2M-1 

■D{Q'■,^^'C = <^'')+ E M^,'md,+,{^,') . (ii) 



j=i 



Here the coefficients djds^/Q"^) are chosen so that the expressions (|31) and ((TU)) are formally identical. The advantage 
of using here the logarithmic derivatives © and the expansions pUj) and ((TTj) .^ as opposed to the expansions (|31) and 
(|^. lies principally in the simple recursion relations for a„'s 



dln/i^ 

whereas for the powers a" the relation is more complicated 

d 



o-nip^) = -Ponan+i , (12) 



d In /i^ 



a{^l'r = -nPoa{p'r+' (l + cM^^^) + C2a{i?f + • • ■ ) , (13) 



the right-hand side here being the consequence of the RGE ([5]). When we use the fact that the full series X'(Q^)mpt 
in Eq. ([TOl) is RScl independent 

'^ ViQ^Upt = , (14) 



(iln/j,2 

we obtain a set of differential equations 

d 



dnip'/Qn^nPodnMfiyQ') (n = I,2,...), (15) 



dhifj? 
whose integration gives (we denote throughout dj{l) = dj and do = 1) 

M^^VQ') = d„ + Y:^ ^j^-^ /3^- In'^ (^g^ j 5-„_. . (16) 

We note that the relations (|16p for dn{^'^ /Q'^), in contrast to those for (i„(/i'^/(5^) in Eqs. (jH])-®, do not involve any 
higher-loop beta coefficients fij (j > 1). Therefore, it is suggestive to compare the situation with the one-loop limit 
of QCD (where Pi = P2 — ■ ■ ■ — ^)- In that limit the perturbative coupling, now denoted as audj,'^), has the one- loop 
RGE running from a given value a{Q^) at the scale Q^ to the scale fj? 



^ Note that the factor in front of the right-hand side is chosen such that Si = a and a„+i = a""*"^ -|- 0(0"+^) for n > 1. Only at one-loop 

level approximation we have a„_|_i = 0"+^, but in general a„_|_i ^ a"+^. 
^ The logarithmic derivatives of the couplin g a nd the expansions of the type IjlOII and l|ll|l were used systematically in Refs. [161 . Il7^ (in 

the context of analytic QCD), and in Ref. [l8| (in the context of pQCD). 



Furthermore, in this case we have a„+i,i^(/^^) = ai£(/i^)"+^, where a„+i^if(/x^) are the logarithmic derivatives of 
QieifJ,^) analogous to Eq. ©. 

Consequently, if we define the (auxiliary) quantity "D^Q") via the following power series: 



2?(Q')pt = a«(M') +^dj(MV(5') aw(M 



2^ + 1 



(18) 



then Eqs. (fT6|) represent the correct /j,^ dependence of the coefficients so as to ensure /i^ independence of the auxil- 
iary quantity T>{Q^). Phrased differently, the auxiliary quantity (|18p is exactly invariant under the combined RScl 
transformations 

5j ~> djifj.'^/Q^) via Eq. HH) , a{Q^) -^ au{t^^) via Eq. ^ . 



(19) 

Note that Eq. (1171) has the form of a homographic transformation. The latter observation leads to an appropriate way 

for treating truncated series, which are in general fj? dependent due to truncation, in particular D((5^;/i^)L (we 
consider truncated series with an even number of terms) . Namely, it is well known in mathematics that the diagonal 
Pade approximants (dPA's), being ratios of two polynomials (Pm, ^Af ), both of order M 



[M/M]{x)^Pm{x)/Rm{X) 
remain dPA's under the homographic transformation 

X *^x = x/{\ + Kx) , 
(where K is an arbitrary constant). This means that 

[M/M]{x)=VM{x)/nM{x) , 



(20) 



(21) 



(22) 



where Pm(x) and TZm(x) are again two polynomials both of order M. More explicitly, if [M/M]f{x) is the dPA of a 



function f{x) whose Taylor expansion around x = exists {f{x) — [M/M]jr{a 



^2M+l 



) , then there exists a function 



F (t^ /) such that [M/M]j{x) = [M/M]f{x). As a consequence, it can be shown that for any function / (with Taylor 
expansion around x = 0) the following identity holds :'^ 



[M/M]fix) = [M/M]j{x) , (23) 

= f{x). In our case of 'D{Q^)pt and its expansion ^TE\\ . we identify: x = a{Q^), 
/3oln(AiVQ'); / = g2 exp(if//3o)2, and V{Q\t = f{x) = f{x). The latter 
identification holds because 2?(Q^)pt — x + Yl°^=i djX^^^ —x + X]?li dj{ii'^/Q'^)x^~^^ . The identity ([23| means that 
dPA's of V{Q^)pt have exact independence of the RScl /i^. Stated differently, when constructing dPA of expansion 



where x = x/{l + Kx) and f{x) 
- auifi^) = x/{l + Kx) [K = 



X 



P^ . it does not matter which value of the RScl n^ we use in ^ 

This fact was noticed by Gardi [19.] , who, as a result, argued that the truncated perturbation series of the form 
^ for physical observables V{Q^) can be well approximated by dPA's because the result is approximately RScl 
independent (i.e., it is exactly RScl- independent when the RGE-running is approximated to be one- loop). Here we 
see that these considerations are valid without approximation for the (RScl-independent) auxiliary quantity I?(Q^) 
which is defined via the power series ^TE\\ . This is related with the fact that the RScl dependence of the coefficients 
dj{ii^ /Q^) as given by Eq. ([TBI) , although involving only /3o and no higher fij coefficients, is exact. On the other hand, 
the RScl dependence of the original coefficients djip^/Q"^) appearing in the power series ^ is more complicated and 
involves (for j > 2) higher-loop beta coefficients /3k {k < j — I), as seen in Eqs. (|6l)-([8]). 

The dPA [M/M] of V(Q'^) has the general form 



[M/M]^{au{p^))--x 



1 + ^12: 



■A 



M-lX 



M-1 



l + Bix 



iM-i 



x=ai{(p2) 



(24) 



We have: f{x) - [M/M]j{x) ~ i^^M+l^ ^^^ jj-^-, _ [m/M]j{x) ~ x^^'+'^ ~ x^^^+i. Therefore, since f(x) = f(x) and [M/M]j{x) = 
[M/M]f{x), we obtain: [M/M\f{x) - [M/M]f{x) ~ x^^f+l. This impUes [M/M]f{x) = [M/M]f(x) (i.e., [M/M]f{x) = [M/M]j(x), 
Eq. II23I I). because the [M/M]{x) Pade's are uniquely determined by the coefficients of their expansion in powers x" for n < 2M. 



We rewrite it by applying a partial fraction decomposition of the fraction on the right-hand side.^ If we denote the 
M zeros of the denominator polynomial (1 + Bix + • • • + Bmx^) by — 1/wj (j = 1, . . . , M), we obtain 



iVl 



M 

(25) 



J = l 

with appropriate "weights" aj (j = 1, . . . , M). Using Eq. ([TTl) gives us finally 

M 

[M/M]^{au{fi^)) = Yl «J a«(Q') ' ^^ere Q^ = ^2 exp(iIj.//3o) , (26) 

i.e., we expressed [Af/Af ]p as a weighted average of one-loop running couplings defined at specific reference momentum 

values (gluon virtualities) Q^ (j = l,...,Af).^ Since, as argued, the expressions (P^ - (|26p are exactly independent 

of the RScl chosen in the original series (IT51) , both the weights cxj and the scales Q? are exactly independent of this 
RScl. 

This observation helps us find an analogous approximant for the true observable T) (or its truncated version 2?P*^1). 
By comparing Eq. (|10p with (IT5|) . we are motivated to define the following approximant: 

M 

er'''(Q^)=E".«(^')' (27) 

i.e., we simply replace in the expression (j26l) the one-loop running coupling aii{Q^j) by the exact (n-loop running, n 

arbitrary) coupling parameter a(Q^). 

The resulting approximant has two important properties: 

1. It is, by sheer construction, exactly RScl invariant (since ctj and Q^ are independent of ^^); 

2. It fulfills the approximation requirement 

V{Q') gl^'/''\Q') = 0(a2M+i) = 0{a'"+') , (28) 

i.e., it reproduces the first 2Af terms of the series PH]) and of the series ([3]). It is relatively straightforward to 
show the latter fact, by expanding the expression (|27|) in terms of logarithmic derivatives (see the Appendix). 

An approximant of the type (j27p was originally introduced in Ref. 14], based on more mathematical considerations. 
It was called "modified Baker-Gammel approximant" and interpreted as a particularly clever resummation procedure 
for the physical observable 'D{Q'^). In Ref. [IJ], also a more formal proof of the properties [1] (RScl invariance) and [2] 
(approximation property) was given. The proof rested on choosing the kernel of the Baker-Gammel approximant to 
be k{z,u) = f{u)/z where z = a{fi^), u = /3o\n{Q'^/ ^'^) and f{u) = a{Q^).^ Within the present paper we constructed 
the same approximant (|27p in a more heuristic and physically motivated manner. 

In Ref. [13 we extended the construction of this approximant so as to be applicable also to the case when an even 
number of coefficients dj (j — I, . . . , 2 A/) are known in the expansion ^, and in Refs. [21| the method was applied 
in pQCD. 

We can interpret the form (j27p as a kind of extension of the previously known scale-setting techniques (principle of 
minimal sensitivity |23| , effective charge method 
Lepage-Mackenzie [26| and its extensions [23, [23, 



2J] and related approaches [25|, and the scale-setting of Brodsky- 
301 1) to several scales. However, in the presented case these scales 
are not fixed by a specific motivated prescription of scale-setting, but are rather based primarily on the successes of 
diagonal Pade approximants in physics and on the additional requirement of refining the approximate (one-loop) RScl 



* In Mathematica [20ll . the command "Apart" achieves this. 

^ In principle, —1/uj's (and thus Q^'s) and 3j's can be sorted into complex conjugate pairs and into real values. In Sec. llVl we apply 

this approach to the massless Adler function for which it turns out that all Q^ and cij are real. 
® For the conventional Baker-Gammel approximants, see for example part II of Ref. [22l | . Exact RScl invariance of such constructions in 

the special case of the aforementioned kernel was apparently first shown in Ref. [14|| . 



invariance of the approximant to the exact RScl invariance. These approximants are global, i.e., they go beyond the 
polynomial form in a, and this is one of the reasons why we expect them to include nonperturbative effects. 

Also interesting to note is the connection of our approximant p7p with Neubert's resummation method [3l| which is 
defined by integration over the momentum flow within the running coupling parameter and the connected momentum 
distribution function w-p 

/•oo 

V^^^\Q%t^ dtwv{t)a{tQ^e'^) . (29) 



Here, C — —5/3 if the "MS" convention for the scale Aqcd is used. When expanding the parameter aitQ^e'-') around 
a{iJ?'), it turns out that this expression represents exactly the leading-/3o part (LB) of the "modified perturbation 
expansion" ((TU)) (cf. Ref. [ij] , and Eq. (PT|) later in the present paper) . We see that our approximant (P7)) is equivalent 
to an approximation of the distribution function wx> (t) in the integrand in ()29p in terms of the weighted sum of delta 
functions 



M 

where the delta peaks are located at tj's such that tjQ^e'^ ~ Q^j {J — ^^ ■ ■ ■ ^ ^)- 

III. APPLICATION TO ANALYTIC QCD MODELS 

In general, the perturbative QCD coupling a{Q^) has a cut in the complex Q^ plane along the negative semiaxis 
up to the positive Landau branching point A^. On the other hand, by the general principles of the local and causal 
quantum field theory [ll, 0, the spacelike observables 'D{Q'^) (such as the Adler function, sum rules, etc.) must be 
analytic functions in the Q^ complex plane with the exception of the cut on the negative semiaxis Q^ G C\(— oo, 0]. 
This analyticity property, however, is not reflected by the a{Q^) which has a cut on a part of the positive axis [0, A^]. 
Therefore, various analytic QCD models have been constructed where the nonanalytic a{Q^) is replaced by an analytic 
Ai{Q^) which has no singularities for Q^ G C\(— oo,0] and at high \Q'^\ ^ A^ (approximately) agrees with a(Q^). 
For details on some of such models we refer to various references: minimal analytic (MA) model 3r5:. 32]: modified 
minimal analytic model [33|; analytic perturbative models |34| : a specific ("close to perturbative") analytic model 
[35!]. Reviews of analytic QCD models are given in Refs. [6|-0|. Calculational techniques applicable to any analytic 
QCD model (the latter being defined via a specification of „4i(Q^) only) are described in Refs. [16, 17, 36]. 

It is natural to ask: how do our approximants (/[*^/*^l fare in such analytic QCD models. As mentioned above, these 
approximants (P7)) choose specific scales which, for low-energy observables, are often close to or inside the (unphysical) 
Landau singularity regime of a{Q^). Therefore, the hope is that our approximants fare much better or even develop 
all their potential in analytic QCD models where they look simply as 

M 

ar''l(Q^an.)=5:5,^i(g,^). (31) 

The other intriguing aspect is that, in any analytic QCD niodeF the analytization of the higher powers a" goes in 
fact via the analytization of the logarithmic derivatives ([9]), cf. Refs. |16l.ll7| 

an+i^^An+i (n = 0,l,2,...) , (32) 

where An+i are the logarithmic derivatives of the analytic coupling Ai 

A.+i(g^)^^5A©, (n = 0,l,2,...), (33) 



We regard the specification of the couphng function Ai{Q^) in the complex Q^ plane as the full specification of an analytic QCD model. 



and not via the naive replacement a" i— J> „4".^ This means that the evaluated observables in analytic QCD have the 
(truncated) "modified analytic" (man) series form analogous to the (truncated) "modified perturbation" (nipt) series 
form in pQCD (fTO jl - ifTTI) 

OO 

2>(Q')ma„ = A^{^i') + J2dJ{^^VQ') A+i(a^') , (34) 

2AI-1 

^(Q';A^')£a? = A^i^,')+ Y. d,{^^^|Q^)A,+,{^?) . (35) 

In view of the presented resummation method (j27l) . this is intriguing, because it shows that the series in logarithmic 
derivatives of the coupling play a central role both in the mentioned resummation method [cf. Eqs. (jlOp . (1181) ] and in 
the evaluation procedure in analytic QCD models [Eqs. (P^ - (|55|) ]. 

The reason for the necessity, in the analytic QCD models, of the evaluation of the observables via Eq. (|55|) originates 
from the fact that the unphysical renormalization scheme (RS) dependence of the truncated series (|35p is^ 



■ = kN{^l'/Q')AN+l{^^') + 0{AN+2){-AN+l), (RS = ln/i^c2;c3;...) , (36) 



OLJyQ ; RS)man _ 7 / 2 //o2\ T (2 



9(RS) 

and from the fact that in analytic QCD models we have the hierarchy Ai{^^) > |^2(/^^)| > lA^di^)] • • • at all complex 
11^. We stress that the expression on the right-hand side of Eq. ([55]) contains only terms ^^(/i^) {j>N + l) and no 

other type of terms. For example, if RS=ln^^, the right-hand side of Eq. (pS)) is exactly —l3QNdff_l{fl'^/Q^)AN+l{^J''^)■ 
If we performed the evaluation by the replacement a" t— > Ai {n > 2), the resulting truncated analytic power series 

2M-1 

P(0^;RS)il'4s = A(m') + E d,{f,'mA^{fiy+' . (37) 

i=i 

would possess in general an increasingly strong RS dependence when the order of the truncation N increases 

^^^^gggj^"^^^ = k^ip'mA^'U^') + 0{Ar') + NPa. , (38) 

where the terms NP^v denote nonperturbative terms {^ (A^//!^)''), which in general become more complicated and 
increase in their value when N increases. The origin of such terms is the difference Ai{fJ,^) — a{fi^) ^ (A^//i^)™ at 

/i2 > A2. 

It is evident that our approximant in analytic QCD, Eq. (|3ip . is RScl invariant (since 5j and Q^ are). Furthermore, 
in complete analogy with the pQCD case, we can show that it fulfills the approximation requirement analogous to 
Eq. (HHD 

V{Q%,,n - Q't"^'^ (Q'; an.) = 0(l2M+i) , (39) 

where the right-hand side has only terms of the form AjiC^'^) {j > 2M -f 1). The relation (|39| . together with the 
aforementioned hierarchy of Aj^& in analytic QCD, gives us additional hope that our approximants pi|) will give us 
values increasingly close to the full value I'((5^)man, Eq- (|M|) . in any chosen analytic QCD model. We will see in the 
next Section, on the example of the Adler function at low momenta (Q^ = 2 GeV ) that this hope is well grounded. 



* The analytic analogs An{Q'^) of powers a(Q^)" are obtained from the relations A-n, = An +5I]m>i km{n)An+m, where the coefficients 

km(n) are obtained from the corresponding pQCD RGE equations (with An i— > a", An+m i— > 'o.n+m)- These relations were presented for 
any analytic QCD model in Refs. [16, 17] in the case of integer n, and in Ref. 1391 for noninteger n = u. The recurrence relations leading 
to the above relations, for integer n and within the context of the minimal analytic (MA) model of Refs. [a-[3. |32|| , were presented in 
Refs. 171. 13711 . Such construction of higher power analogs An, not as powers of A\ but rather as linear (in A\) operations on A\, reflects 
a very desirable functional feature: their compatibility with linear integral transformations (such as Fourier or Laplace) [331 ■ On the 
other hand, in linear tranformations, the image of a power is in general not the power of the image. 
^ The relation I I36II can be obtained in complete analogy with the perturbative QCD, under the correspondence II32I I. 



TABLE I: The weight coefficients a_, and the scale ratios Q^j/Q^ for our RScl- invariant approximants, Eqs. (|27p and l|3ip . for 
various order indices (M = 1,2,3,4), in the case of leading- /3o massless Adler function. 



M 



(5i; QllQ'') 



(52,Qi/Q^) 



(5«,Q3/Q') 



(54,Ql/Q' 



M = l 


(1; 0.5001) 


- 


- 


M = 2 


(0.6948; 0.1711) 


(0.3052; 5.771) 


- 


M = 3 


(0.3579; 0.07969) 


(0.6011; 1.0534) 


(0.0410; 85.77) 


M = 4 


(0.1376; 0.03803) 


(0.6821; 0.3862) 


(0.1767; 17.16) 



(0.0037; 1518.; 



IV. NUMERICAL CHECKS OF THE QUALITY OF THE APPROXIMANTS 



In this Section we will investigate how our approximants (I5T|) [and (P7| ] work when applied to a spacelike QCD 
observable whose perturbation series is known to a sufficiently high order. Specifically, we will consider the massless 
Adler function 'D{Q^) at low Q^ (Q^ — 2 GeV ) and perform numerical evaluations of our approximants both in 
perturbative QCD (pQCD) and in three different analytic QCD (anQCD) models, namely: 

• Minimal Analytic (MA) model of Refs. d-d; 

• the approximately perturbative anQCD model of Ref. [3a] (CCEM); 

• the perturbative anQCD model type "EE" (whose beta function involves exponential functions) in two variants, 
of Ref. [13. 

The characteristics of these different models will be specified in more detail later in this Section. Beforehand, we 
sketch the general procedure: we will consider first the leading-/3o (LB) resummation part of V whose expression in 
pQCD is 



dt 



V^^^\Q\, = / ^ Fj,{t)a{tQ'e^) 



a(Q') + di,i/3oS2(Q') 



dn,nPo CLn+l{Q ) + 



(40) 
(41) 



Here, F-pit) = WT>(t)t is the characteristic function of the Adler function, whose explicit expression was obtained in 



Ref. 131| on the basis of the leading-/3o expansion coefficients d 



(LB) 



_^LB) 



— dn.npQ obtained from the leading-/3o 



Borel transform of Refs. [33, EO] at RScl fi^ = Q^ in the "V" scale convention. ^° The coefficient d„ ' represents 



simultaneously the leading-/3o part of d„ and of dn once these two coefficients are organized in series in powers of nf 

and thus of /3o; dn is RSch independent but RScl dependent (see also Eq. ([T6l) : for details, see Ref. 17|). 

The evaluations will be performed in the simplest renormalization scheme C2 = C3 = • • • = in various QCD 
models (pQCD and anQCD's, except the anQCD model "EE"). This is convenient because the expressions are 
then simple and explicitly related with the Lambert function [l2|, |41| . As the point of reference we take the value 
a(M|, MS) = O.HQ/tt. This then corresponds to the value a{p,f^]nf = 3; C2 = C3 = • • • = 0) « 0.2215/7r at the 
"initial" chosen scale /ijn — 3toc — 3.81 GeV. 

We will assume that Uf — 3 in our calculations. At Q^ — 2GeV^ we obtain a(2GeV'^) — 0.3479/7r. 

The practical evaluations can be performed by choosing any value of RScl fj,^, e.g. fj? — Q^. In the leading- /3o 
case the choice fi^ = Q^ means using the coefficients (i„^„ = (i„^„(/i^/(5^ = 1) in the expansion (|411) . Nonetheless, 
as shown, the use of different RScl fj,'^ ^ Q^ gives us identical results, as can be checked numerically as well. We 



note that by choosing /i = Q , the coefficients dj = dj(l), dj = dj{l) are Q -independent, 
coefficients cij and parameters Uj in Eqs. (I?7l) and (I5T]) are Q^-independent (when fj.'^ = Q^] 



Therefore, the weight 
, and thus the ratio of 

scales Q^j/Q'^ = exp{uj/(3o) [see Eq. p6)) . with /i^ — Q^] will be Q^-independent (and, of course, ^^-independent). In 
TableUwe give the values of weights cij and scale ratios Q'^/Q'^ for various indices M of our approximants. We can see 

from the Table that the scale ratios Q^/Q^ get increasingly spread out when the order index M increases. However, 
for those ratios which are much smaller or much larger than unity, the corresponding weight factors are small. 



Ref. 13111 uses the notation D{t) = 4Fx>{t)/t. Note that we use throughout the "MS" convention for the scale A, i.e., C = C = —5/3. 
Large-/3o calculations _are usuallj;^ performed with "V" scale convention, Le., C = 0. The relations between the two, at a given RScl fj? 



(e.g., n^ = Q^), are: d„,„(A) : 



.(Av) + Efc=i (n!/(fc!('^ - fc)!)(-C)'=rf„ 



-fc(Av) + (-C)". 



TABLE II: The results of the one-loop approach: diagonal Fade approximants (dPA) ([M/M]) with increasing index M, at 
Q^ — 2 GeV^, for the leading-/3o massless Adler function 'D(Q^). For comparison, the result of the Principal Value of integration 
(with the estimated IR renormalon ambiguity) is included. In addition, the approximants (|27|) in the case of full pQCD running 
a{tQ^e ) are included, and the corresponding Principal Value. In the parentheses, the corresponding results of the truncated 
series (|4H) are given (with RScl /i^ = Q^). See the text for details. 



Case M = 1 M = 2 M = 3 M = 4 M = 5 M = 6 



PV 



1-loop 0.134(0.130) 0.161(0.155) 0.175(0.164) 0.194(0.160) -0.497(0.080) 0.156(-0.714) 

full 0.140(0.134) 0.200(0.174) 0.532(0.198) 0.095-0.051i(0.107) 0.162-0.009i (-1.79) 0.250-O.OOli (-39.i 



0.178 ±0.020 
0.174 ±0.020 



The authors of Ref. [44] applied the diagonal Fade approximants to the (auxiliary) power series quantity P^'"^) (Q^)pt 
(at Q^ — 2 GeV^) obtained from the series (|4T|) by the replacement a„+i i-^> a"+^ (the approximation of one-loop 
RGE running), and compared with the result of the integration (j40l) obtained by assuming one- loop RGE running 

of aitQ^e'-'); the integral is ambiguous in the integration at low t (IR regime) due to the Landau singularity, so they 
chose the Principal Value for the integration. 

The results of this type of (one- loop) evaluation are given in Table HIl for the case Q^ ~ 2 GeV^. We fix the 
one- loop running coupling au{Q ^) so that it agrees with the aforementioned full a at Q^ = 2 GeV^: au{Q'^) — 
o-iQ"^) = 0.3479/7r. In addition, we include in the Table the corresponding results with the full pQCD evaluation in 
the C2 = C3 = • • • renormalization scheme ("two-loop") which uses in the integral PO]) the full pQCD a{tQ^e'^), and 
our approximants (l?7l) . We can see that the dPA's (in the one- loop case) and our approximants ((Tf)) oscillate rather 
erratically around the corresponding Principal Value. This has to do with the fact that, at higher order index Af 
(Af > 3) the scales Q^ come rather close to the Landau singularity of the running perturbative coupling. In fact, 

the approximants become even complex in the full case once at least one of the scales Q^ hits the unphysical cut 
(0,Al ) (where: Al ~ 0.150 GeV^, i.e., Al. ~ .388 GeV), since a{Qj) becomes complex. In the one-loop case, we 
have a simple Landau pole instead of the cut (with A^ ~ 0.036 GeV^, i.e., Al. ~ 0.190 GeV), so the approximants 
would remain real even when one of the scales were below the Landau pole. In the parentheses, the results of the 
corresponding truncated series are given - for the one-loop case the truncated version V{Q^; /^^)pt of the expansion 

([T8|. and in the fuU (loop) case the truncated version 2?((5^; /i^)l„pt Eq- pH . both with RScl ^^ = Q^. We see 
that these truncated series behave in general worse than the resummed versions, and show for larger Ad asymptotic 
divergence (in the one-loop case for M > 5, and in he full loop case for Af > 4). 

There are several analytic QCD models (for AiiQ"^)) in the literature. The most used one is the model of Shirkov, 
Solovtsov and Milton [S-IH, which keeps for the cut of Ai{Q'^) on the negative Q^ axis the discontinuity function of 
the pQCD coupling a{Q^), and the unphysical pQCD cut on the positive axis is eliminated 

A^HQ')-- da^, (42) 

where p^' (cr) = Ima((5 ^ = — c — *e)- This represents, in a sense, the minimal changes (in the cut) with respect to 
pQCD. Therefore, we call this model the Minimal Analytic (MA).^^ The only adjustable parameter there is the scale 
A (in the "MS" scale convention). In order to reproduce QCD phenomenology at high energies, the value of this scale 
at n/ = 5 in MA is about 260 MeV, which corresponds at ny = 3 to the value of A w 415 MeV [9|. We will use this 
value in MA, and will use there also the RSch C2 — c^ = ■ ■ ■ = 0. 

Another analytic QCD model is described in Ref. [35] (CCEM). It differs from MA in the sense that the discontinuity 
function pi(cr) — Ini^i(— cr — ie) differs from the pQCD discontinuity function at low cr < 1 GeV where it is replaced 
by a delta function. The spacelike coupling Ai is then 

Ai{Q') = -^ + - r ds'-£^ , (43) 



^'^ In the literature, it is usually called Analjrtic Perturbation Theory (APT), and it then involves a specific construction of the analytic 
analogs of higher powers a". The construction can be applied only in MA, and it is in such a case equivalent to the construction 



presented in Refs. llft llTll (the latter construction being applicable to any anQCD). 
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where u = Q^/A^^, s = a/K^, r'^^^Xs) = p^^^\a) (in the RSch ca = C3 = ■ • • = 0), and Aw « 0.487 GeV 
is the scale appearing in the Lambert function W:fi{z±). The scale Aw was fixed basically by the requirement 
that the high energy QCD phenomenology be reproduced. The (dimensionless) iree parameters (/f , si = Aff /A^, 
So = Mq/A^^) are fixed in the model in such a way that at high Q^ it merges with the pQCD coupling to a high 
degree of accuracy [AiiQ"^) — a{Q^) ~ (A^/Q^)^] and that, simultaneously, it reproduces the measured value of the 
semihadronic (massless and strangeless) tau decay ratio^^ rr{AS = 0,mq = O)oxp. = 0.203 ± 0.004. We note that 

in MA we have Ai (Q^) — cl{Q^) ^ (A^/Q^), i.e., at high energies this difference is not quite negligible, and the 
predicted value of ^^(AS' = O^iriq — 0) is about 0.14. 

Yet another analytic QCD model which we will use is the so called EE model of Ref . 34] , which is in fact a fully 
perturbative analytic QCD model [the /3(a) function is analytic function of Ai{Q'^) = a{Q^) at a = 0].^"^ The beta 
function has the Ansatz 

/3(a) = -/3oa2(l - y)/(r)|y,,/a„ , (44) 

where ao = a{Q'^ = 0) is a finite value (infrared fixed point), f{Y) is analytic at F = 0, and we require analyticity of 
a(Q^) at Q^ = 0, which turns out to give the condition ao/3o/(l) = 1. The expansion of /3(a) in powers of a also has to 
reproduce the first two universal coefficients /3o and /3i , cf. Eq. ([5]) . There are at least two variants of the mentioned 
"EE" model. In the first variant ("EEvl") the function f{Y) in the beta function is a combination of (rescaled and 
translated) functions (e^' - \)/Y and Y/{e^ - 1) {e^ - l)/Y and Y/{e^ - 1): 

^^ ' [fci(y-yi)] (exp[-fc2(i"->'2)]-i) ^1' 1' 2- 2^, V ^ 

where the constant JC ensures the required normalization f(Y = 0) = 1. In this variant we have, at first, five real 
parameters: ao = a{Q^ = 0) and Yj, kj {j = 1, 2). Two parameters (Y2 and ao) are eliminated by the aforementioned 
conditions: ao/3o/(l) = 1 and the reproduction of the universal /3i coefficient. The other three parameters are 
approximately fixed by the condition of analyticity of a{Q^) and the requirement of obtaining as high a value of 
rriAS = 0,TOg = 0) as possible (it is always too low in comparison to the experimental value 0.203 ± 0.004). The 
obtained values are: Yi = 0.1, ki = 10.0, fca = 11.0. This results in ao = 0.236 and the highest possible value 
rr{AS = 0,TOq = 0) « 0.15. This latter value is still clearly too low. 

The second version ("EEv2") has the function f{Y) in the beta function modified, in comparison to EEvl, by a 
factor /fact 

EEv2 : /eev2(>") = hE.i{Y)fi,ct{Y) , (46) 

with /facUr) - (1 i\B +]P)Y^) (1 « ^ « ^) ■ (47) 

This factor has the values of K and B adjusted so that the expansion of the evaluation of r^iAS = 0, iriq = 0), by 
the inclusion of the leading-/3o (LB) contribution and of the first three beyond-the-leading-/3o (bLB) contributions, 
gives the correct r^ value: rr{AS = 0,m, = 0) = 0.203 (=> i? = 1000 and K = 5.4). The factor /fact(5^) does not 
destroy the analyticity of a{Q^), and it does not change substantially the values of a{Q^) since it is close to the value 
one for most F's. However, the price that we pay is high nonetheless: the coefficients Cj = Pj/Pq of the expansion 
of the modified beta function are extremely high for j > 4 {cj > 10^ for j > 4), implying strong divergence of any 
evaluation series of observables (including r-r) when bLB terms of ^ a" with n > 5 are included. The factor /fact(Y) 
introduces singularities of /3(a) at rather small values of \a\. 

For more details on the models CCEM (with so = 3.858) and EEvl and EEv2, we refer to Refs. @ and (H], 
respectively. 

The results of our approximants (I5T|) in these analytic QCD models, for the lcading-/3o part of the Adler function 
V{Q^) at Q^ = 2 GeV^, are presented in Table Hill For comparison, the exact integrated values 

-D^fKQ') = f jF^{t)A,{tQ'e^) (48) 



^^ We use the variant of the model with the value of sq = 3.858, which reproduces the measured value of Vr when the leading-/3o 

resummation and the inclusion of the known beyond-the-leading-/3o terms is performed in the evaluation of r-r- 
^^ Our general construction of An{Q'^) gives in such models: A.n = ^", as it should be. 
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TABLE III: Evaluations of the leading-/3o massless Adler function 7)^^^' (Q^) in various analytic QCD models, using our RScl 
invariant approximants (|3H) . with increasing index M, at Q^ — 2 GeV^. For comparison, the exact result of the integration 
(|48|l is included. In parentheses in the Table, the values of the corresponding truncated series V'^ \Q^)[nan are given (with 
RScl fi^ = Q^). See the text for details. 



model M — 1 



M 



M 



M 



M 



M = 6 



M = 7 



MA 0.1167(0.1147) 0.1222(0.1214) 0.1217(0.1208) 0.1217(0.1205) 0.1217 (0.1211) 0.1217(0.1209) 0.1217(0.1174) 

CCEM 0.1371(0.1321) 0.1649(0.1640) 0.1650(0.1733) 0.1617(0.1788) 0.1624(-0.0048) 0.1632(-0.0407) 0.1626(11.70) 

EEvl 0.1062(0.1047) 0.1141(0.1144) 0.1136(0.1146) 0.1131(0.1138) 0.1132(0.1063) 0.1133(0.0842) 0.1133(12.42) 

EEv2 0.0965(0.0952) 0.1036(0.1035) 0.1035(0.1041) 0.1032(0.1041) 0.1032(0.1018) 0.1032(0.0840) 0.1032(-0.4615) 



exact 



0.1217 
0.1627 
0.1133 
0.1032 



TABLE IV: Evaluations of the full massless Adler function in various analytic QCD models, using our RScl invariant approx- 
imants pi|) for M — 1,2, at Q^ = 2 GeV^. For comparison, two other evaluations (LB+bLB; and tman: truncated modified 
analytic series) are included. See the text for details. 



model 



M 



M ■ 



LB+bLB 



tman 



MA 
CCEM 
EEvl 
EEv2 



0.1175 
0.1389 
0.1070 
0.0972 



0.1196 
0.1535 
0.1164 
0.1390 



0.1191 
0.1528 
0.1183 

0.1584 



0.1199 
0.1541 
0.1195 
0.1587 



are also given in the Table. Note that the leading- /3o integration, Eq. (US)), has now no ambiguities since no Landau 
singularities exist, in contrast to the pQCD case (j40p . Incidentally, the expansion of Eq. (H5|) is completely analogous 
to the pQCD expansion (HIT) 



P^^^)(Q')man = ^i(Q') + rfl,l/3o^2(Q') 



dn^nPoAn+liQ ) 



(49) 



In parentheses, we give the results of the corresponding truncated version of the series (H^l) . i.e., I?(^^) ((9^)man , with 
fi^ = Q^, for each M. We see in the Table that our approximants converge systematically and fast to the exact values 
when the order index M increases. The truncated series, on the other hand, have divergent behavior which, though, 
starts manifesting itself at large M's {M > 7 in the MA case; Af > 5 in the CCEM and EE cases) since these are 
analytic QCD models. Despite this divergence, the aforementioned hierarchy of the couplings |^fe((9^)| > |-4fc+i(Q^)| 
in general turns out to be true for all relevant indices k in the Table (fc = 1, . . . , 13), at Q^ = 2 GeV^. 

The first three coefficients dj {j — 1, 2, 3) are now exactly known for the Adler function p3l - |4a |. Therefore, we can 
construct our approximants (PT|) for the order indices M — 1 and Af = 2 on the basis of these exact four coefficients. 
The results of this calculation, for the three analytic QCD models, are presented in Table IIVI For comparison, we 
also include the results (tman) of the truncated modified analytic series ((35|) . with jjp' — Q^ , and the more refined 
"LB-f bLB" evaluation which takes into account the leading-/3o resummation contribution pS)) and the three additional 
known terms (bLB: beyond-the-leading-/3o) 



^(bLB)(g2 



3 

E 



(d„-d„,„/3?)A+i(Q'). 



(50) 



We can see that our approximants pi|) . with index M = 2, represent a competitive evaluation of the observable, 
especially when comparing with the (partially) resummed results LB+bLB and the truncated (modified) analytic 
series (tman). 

The results of our method with M = 2, in the MA and CCEM cases, deviate from the LB+bLB results less than 
the tman results deviate. Since the analytic models MA and CCEM are in "tame" RSch's [i.e., the ones where the 
RSch parameters Cj (j > 2) are very small, in fact, zero], we can expect that both the LB+bLB and tman approaches 
give good estimates of the true value in the model, and that LB+bLB is probably a better approach since it uses 
significantly more input information than tman. However, we recall that our M — 2 approximants use as little 
input information as the truncated (tman) approach, i.e., the first three dj's, and yet Table HVl indicates that our 
approximants with M = 2 are competitive with the LB+bLB approach in the MA and CCEM models. 

On the other hand, the RSch coefficients Cj are increasing fast in the models EEvl and dramatically fast in EEv2. 
In that case, the coefficients dj and {dj — dj^P^), which depend on Cj via an additive term —Cj/{j — 1) (if j > 2), 
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increase very fast when j increases, so that tman and LB+bLB approaches become uncertain. ^^ We notice that in 
the case of EEv2, our approximant (for M = 2) is essentially different from the LB+bLB and from the tman result. 
The tman series (|35|) and the truncated bLB series (f50|) become in that case very divergent once we include the terms 
An+i with n > 4 (cf. Ref. [3J| for further details on the divergence of the coefficients d„ in this case). In that case, 
our approximants, for M = 2, are probably comparatively the most reliable estimate of the true result in the EEv2 
model. 

V. CONCLUSIONS 

We tested in various analytic QCD models an earlier developed |14l . Il5| RScl invariant resummation method, by 
applying it to the evaluation of the massless Adler function V{Q^) at low energy {Q^ — 2 GeV^). The method 
is global, i.e., nonpolynomial in the (analytic) coupling parameter. It is related with the method of diagonal Fade 
approximants (dPA's), representing an extension of the dPA method by achieving exact RScl independence. The 
method, applied to spacelike observables, results in a linear combination of coupling parameters at several spacelike 
momentum scales (each of them RScl invariant), and thus represents an extension of the well-known scale-setting 
techniques of Stevenson [23| , Grunberg [2J] , and Brodsky-Lepage-Mackenzie [2a] . For observables with low scale Q^ 
of the process, the method when applied within the perturbative QCD is not very efficient in practice. The reason 
for this is that the perturbative QCD coupling a{Q^) has unphysical (Landau) singularities at low positive Q^, and 
some of the scales of our approximant turn out to be close or even within this singularity sector. On the other hand, 
the method turns out to be very efficient in analytic QCD models, because the analytic coupling Ai{Q'^) has no 
unphysical singularities. In the case of the leading-/3o part of the Adler function, the results of the method converge 
very fast to the exact result within each analytic QCD model. Furthermore, when the method is applied to the 
truncated (analytic) series of the entire Adler function, whose first three coefficients beyond the leading order are 
known exactly, the result of the method becomes competitive with the result of the sum of the (exact) leading-/3o (LB) 
contribution and the truncated beyond-the-leading-/3o (bLB) analytic series, although the latter method (LB-|-bLB) 
uses significanly more input information than our method. We conclude that our method is at the moment probably 
the best method, in the analytic QCD frameworks, for the evaluation of spacelike observables when the evaluation is 
based on the known part of the truncated integer power perturbation series of the observable. The method can be 
used also for the evaluation of timelike observables (such as the cross section of e~^e~ scattering into hadrons, and 
semihadronic r decay ratio Tt) when the latter are expressed as contour integrals involving spacelike observables. 
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Appendix A: The approximation requirement 

Here we show that the approximation requirement, Eq. (j28p . is fulfilled by our approximant (|27p. Taylor-expanding 
a((5^)'s in the approximant around ln(/i^), by using the definitions ©, we obtain 

M M 



j = i i=i fc=0 

oo AI , oc AI 

= J2ak+i{fi')J2a, (-/3oln(Q2/^2)^ ' ^ Y,ak+i{^^')Y.S,{~u,f . (A2) 



fc=0 j=l k=0 j = l 



14 The Cj = Pj/Po coefficients in EEv2 are: -106. 8(j = 2); 326. 7(j = 3); 1.72 • 10^{i = 4); 3.08 ■ lO'^O' = 5), etc. 
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In the last equation we used the fact that Uj = /3o ln(Q|/^2), see. Eqs. ((251) -dM])- However, Eqs. (I24l)-(l26l) and p8)) 
teU us that 

ViQ%t - [M/M]^ {au{^i'')) = O {au{^l^f^'+') . (A3) 

This imphes that the expansion of [M/M]g(x) in powers of x = ai£(/i^) reproduces"'^^ the coefficients at powers of x" 
for n = 1, . . . , 2M in the expansion of V{Q'^), Eq. ^ 



M 

'^a^i^Uj)'' ^dk{n^/Q^) ior k = 0,l,...,2M -1 . (A4) 

Note that do{fi'^/Q'^) = 1. Inserting the indentities (jA4l) into Eq. (jA2p . we obtain 

2J\/-1 

a|f /'"(Q') = E afc+i(A*')4(A*VQ') + 0{a2M+i) . (A5) 

fc=0 

This, in combination with the expansion ([TOl) of the observable T>{Q'^) in afc+i(/i^), gives us immediately 

2?(g')mpt - ei''/'"(Q') = 0(52M+i) = 0(a2*^+i) , (A6) 

i.e., the approximation identity (j28l) . The same proof can be repeated in analytic QCD models (except for the 
notational change Ofc+i *-^ •A.k+i), i.e., the approximation identity p9l) is also valid. 
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